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ABSTRACT 

We explore the nature of Infrared Excess sources (IRX), which are proposed as candidates for 
luminous [L x (2 - lOkeV) > 10 43 ergs" 1 ] Compton Thick (N H > 2 x 10 24 cm- 2 ) QSOs 
at z w 2. Lower redshift, z « 1, analogues of the distant IRX population are identified by 
firstly redshifting to z = 2 the SEDs of all sources with secure spectroscopic redshifts in the 
AEGIS (6488) and the GOODS-North (1784) surveys and then selecting those that qualify 
as IRX sources at that redshift. A total of 19 galaxies are selected. The mean redshift of the 
sample is z « 1. We do not find strong evidence for Compton Thick QSOs in the sample. 
For 9 sources with X-ray counterparts, the X-ray spectra are consistent with Compton Thin 
AGN. Only 3 of them show tentative evidence for Compton Thick obscuration. The SEDs of 
the X-ray undetected population are consistent with starburst activity. There is no evidence 
for a hot dust component at the mid-infrared associated with AGN heated dust. If the X-ray 
undetected sources host AGN, an upper limit of Lx(2 — 10 keV) = 10 43 erg s -1 is estimated 
for their intrinsic luminosity. We propose that a large fraction of the z ps 2 IRX population 
are not Compton Thick QSOs but low luminosity [Lx (2 — 10 kcV) < 10 43 erg s _1 ], possibly 
Compton Thin, AGN or dusty starbursts. It is shown that the decomposition of the AGN and 
starburst contribution to the mid-IR is essential for interpreting the nature of this population, 
as star-formation may dominate this wavelength regime. 
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1 INTRODUCTION 

The composition of the diffuse X-ray background remains a prob- 
lem for high energy astrophysics. X-ray imaging surveys with the 
Chandra and XMM-Newton observatories have resolved between 
80 to 100% of the XRB into discrete sources below about 6keV 



(e.g. Worsley et al. 2005 , Hickox & Markevitch 2006; Georgakakis 



|et al.|2008 1. The vast majority of these sources are Compt on thin 



AGN (N H < 10 cm -2 ) at a mean redshift z « 1 (e.g. Barger 



et al. 2005) |Akylas et al.||20 06). At higher energies however, be 
tween 20 and 30keV, where the bulk of the XRB energy is emitted 
(Marshall et al. 1980 1, only a small fraction of its total intensity has 
been resolved into discreet sources ( Sazonov et al. 2007). As a re- 
sult the nature of the populations that make up the XRB at these en- 
ergies is still not well known. Population synthesis models use our 
knowledge on the properties of the X-ray sources below ~ lOkeV 
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to make predictions on the nature of the X-ray populations close to 
the peak of the XRB (e.g. |Gilli et al.|2007) . These models indicate 
that Compton thin AGN alone cannot account for the shape of the 
XRB spectrum at ~ 20 — 30 ke V. An additional population of heav- 
ily obscured, Compton thick (Nh > 10 24 cm -2 ) AGN is postulated 
to reconcile the discrepancy (e.g. |Gilli et al.|2007) . The required 
number density of such sources is however, under debate ( |Treister| 
|et al.|2009a]|Draper & Ballantyne|2009| >. Unfortunately the identi- 
fication of the heavily obscured AGN population predicted by the 
models is far from trivial. The X-ray emission of these sources be- 
low about lOkeV is suppressed by photoelectric absorption and as 
a result most of them are expected to lie well below the sensitivity 
limits of the deepest current X-ray observations. Although a hand- 
ful of Compton thick AGN candidates have been identified in deep 
X-ray surveys ( [Tozzi et al.|2 006 ; Georgantopoulos et al. 2009jl, the 
bulk of this population remains to be discovered. 

Selection at the mid-IR (3 — 30 /im) is proposed as a powerful 
tool for finding heavily obscured X-ray faint AGN. The UV/optical 
photons emitted by the central engine are absorbed by the gas and 
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dust clouds and appear as thermal radiation with a broad bump in 
the mid-IR (vj u units; e.g. |Elvis et al.|1994{|Prieto et al.|2009f . Di- 
verse selection methods have been developed to identify this AGN 
spectral signature in the mid-infra red. | Lacy et al.| ( |2004) , |Stern et al.| 
( 2005 ) and Hatziminaoglou et al. ]( |2005[ propose simple colour cuts 
based on the mid-IR colours of luminous high redshift QSOs and/or 
type-2 Seyferts. Polletta et al. (20061 and Rowan-Robinson et al. 
( 2009 1 fit templates to the broad-band photometry from UV to the 
far-IR. |Alonso-Herrero et al.| < [2006] > and |Donley et al.| < [2007] > se- 
lect sources with power-law Spectral Energy Distributions (SEDs) 
in the mid-IR. The methods above have merits and shortcomings. 
Selection by mid-IR colour for example, is simple but suffers con- 
tamination from star-forming galaxies if applied to the deepest mid- 
IR samples currently available ( Georgantopoulos et al. 2008; Don- 
ley et al. 2008). Template fits are powerful but require high quality 
photometry over a wide wavelength baseline for meaningful con- 
straints. Power-law SEDs provide the most clean samples of in- 
frared selected AGN, but are sensitive only to the most luminous 
and hence, rare sources. 

A much promising method that is believed to be efficient in 
identifying heavily obscured, possibly Compton Thick, AGN is 
based on the selection of sources that are faint at optical and bright 
at mid-IR wavelengths. In its simplest version this method applies a 
cut in the 24 /im over i?-band flux density ratio, /24 /Jm / /r > 1000 
to identify Dust Obscured Galaxies (DOGs; |Dey et al.|2008) at a 
mean redshift z ~ 2. In addition to the limit above IFiore et al.l 
(2008 . 20091 also use the colour cut R—K > 4.5 to select Infrared- 
Excess sources (IRXs). This colour selection is motivated by the 
observational result that redder sources include a higher fraction of 
obscured AGN (e.g. |Brusa et al.|[2005| l. A different approach has 
been adopted by Daddi et al. ( 2007 ). They select BzK sources that 
show excess mid-IR emission relative to that expected based on the 
rates of star formation measured from shorter wavelengths. Despite 
differences in the adopted criteria, all the studies above identify 
a population of galaxies with similar properties in terms of aver- 
age redshift (z ~ 2), mean X-ray properties and mid-IR luminosi- 
ties. These sourc es are believed to be massive galaxies (M s 
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experience intense bursts of star-formation and rapid supermassive 
black hole growth, possibly triggered by mergers (e.g. |Bussmarin| 
et al. 2009b; Narayanan et al. 2009). The apparently brighter sub- 
set of this population (S24 > 300 /J.Jy) are proposed as descen- 
dants of Submillimeter Galaxies on their way to becoming unob- 
scured QSOs (e.g. Bussmann et al. 2009b a , Narayanan et al. 2009 ), 
which will eventually evolve into present-day 4 L* galaxies ( |Dey| 
|et al.|2008||Brodwin et al.|2008^ . Although these sources are un- 
doubtedly important for understanding the co-evolution of galaxies 
and SBHs at high redshift, in this paper we focus on their signifi- 
cance to XRB studies by putting into test claims that they include 
a large fraction of Compton Thick AGN. The two key properties 
of these sources (hereafter referred to as IRX sources), which are 
interpreted as evidence for Compton Thick AGN are (i) their hard 
mean X-ray spectrum and (ii) their faintness at X-ray wavelengths 
relative to the mid-IR. The first point is demonstrated in Figure [T] 
which shows the that the mean hardness ratio of the IRX sources is 
consistent with that of the local Compton Thick AGN NGC 1068 
l |Matt et al7||1999| >. With respect to the second point, the major- 
ity of the IRX sources (80 per cent; Georgantop oulos et al.|2008} 
are not detected in the deepest X-ray surveys available and their 
mean X-ray properties can only be studied through stacking analy- 
sis. Figure [2] shows the average X-ray to mid-IR luminosity ratio, 
Lx (2 — 10 keV) / z/L„5.8/im, of the X-ray undetected IRX popu- 



lation. Relative to local AGN from the sample of |Lutz et al.H2004| > 
these sources appear underluminous in the mid-IR by 2-3 dex. As- 
suming a narrow range of Lx (2 — 10 keV) / ^L„5.8/im for typical 
AGN (e.g. Lutz et al. 2004), the position of the IRX population in 
Figure|2]is consistent with a Compton Thick obscuring screen that 
suppresses the observed X-ray emission relative to the mid-IR. Un- 
der the assumption that the mid-IR luminosity is a good proxy of 
the AGN power, one can estimate intrinsic 2-1 0keV luminosities 
for the IRX sources in excess of 10 43 ergs _1 1 Daddi et al.|2007} 
|Fiore et al.|2008| [20091 |Treister et al.|2009b| >. AGN with intrinsic 
luminosities above this limit are hereafter referred to as QSOs. IRX 
sources are therefore, excellent candidates for the Compton Thick 
QSOs needed by population synthesis models to explain the XRB 
spectrum at 20-30 keV. Simulations indeed show that the bulk of 
these sources (80 per cent, Fiore et al. 2008; 95 per cent, Fiore et 
al. 2009) are likely to be Compton Thick QSOs, while their es- 
timated space densities are consistent with the predictions of the 
XRB models. 

Although there is no doubt that the IRX population includes 
a fraction of heavily obscured AGN ( Georgantopoulos et al. 2009 ), 
the evidence for Compton Thick sources is far from conclusive. The 
two key properties of the IRX sources, hardness ratio and X-ray 
to mid-IR luminosity ratio, are also consistent with Compton thin 
AGN of lower luminosity, hereafter defined as Lx (2 — 10 keV) < 
10 43 ergs~\ Fi gure[T|for example, shows that the mean hardness 
ratio of this population could be due to moderate obscuring col- 
umn densities of few times 10 23 cm -2 . Moreover, Figure[2]demon- 
strates that lower luminosity Compton Thin AGN (Terashima et al. 
2002) have very low X-ray to mid-IR luminosity ratios, i.e. sim- 
ilar to those observed for IRX sources. This is not surprising as 
both AGN and star-formation contribute to the mid-IR, thereby re- 
sulting to a broad Lx{2 — 10 keV)/i>Lv5.8[im distribution. One 
has to isolate the AGN component at the mid-IR part of the SED 
to get a tight correlation with X-ray luminosity (e.g. Prie to et aT] 
|2009[ >. |Lutz et al.| ( |2004[ l for example, accomplished that by select- 
ing only sources where the AGN dominates in the mid-IR, i.e. those 
without a significant starburst component relative to the AGN. For 
the high redshift IRX sources it is not clear what fraction of the 
mid-IR emission is from AGN heated dust. Only under the strong 
assumption that the bulk the mid-IR luminosity is associated with 
reprocessed radiation from accretion on the central SBH, can one 
infer that these sources are Compton Thick QSOs. A number of 
studies on the IRX sources have already suggested that this popula- 
tion may include a large fraction of lower luminosity Compton thin 
AGN or even dusty starbursts (e.g. (Georgantopoulos et al. 1(2008} 



|Donley et al.|2008[|Pope~et al.|200 8 Murph y~eTaT.|2009| >. 

Elucidating the nature of IRX sources, Compton Thick QSOs, 
Compton Thin lower luminosity AGN or starbursts, requires very 
deep X-ray data and/or AGN/starburst decomposition at the mid- 
IR. Unfortunately, IRX samples at z ~ 2 are too faint at almost 
any wavelength, expect the mid-IR, for such a detailed study. In 
this paper we address this issue by selecting lower redshift (z ~ 1) 
analogues of the distant (z ~ 2) IRX population. The advantage of 
this approach is that the selected sources are apparently bright and 
therefore their SEDs can be constrained over a wide wavelength 
baseline, from the far-UV to the far-IR. The sources are selected 
in fields with some of the deepest X-ray observations available 
(Chandra Deep Field North, All wavelength Extended Groth strip 
International Survey: AEGIS), thereby allowing study of their X- 
ray properties via spectral analysis. An additional advantage of our 
strategy is that there is only about 2.5 Gyr difference in the age of 
the Universe between z — 1 and z — 2 (Ho = 70 km s~ Mpc - , 
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Figure 1. Hardness ratio against redshift. The hardness ratio is defined as 
HR=(H-S)/(H+S), where H and S are the counts in the 1.5-6 and 0.3-1.5 keV 
energy intervals respectively. The (blue) dot is the mean hardness ratio of X- 
ray undetected sources in the COSMOS survey (Fiore et al. 2009) estimated 
using stacking analysis. The vertical errorbar is the lcr uncertainty of the 
HR, while the horizontal errorbar shows the redshift range of the Fiore et 
al. (2009) IRX population. The dotted line (black) shows the expected HR 
of the local Compton Thick AGN NGC 1068 Matt et al. (1999). The (red) 
dashed curve corresponds to the HR of a power-law X-ray spectrum with 
photon index T = 1.9, which is absorbed by a column density of Njj = 
1.5 X 10 23 cm -2 . The continuous (black) line corresponds to the HR of 
an unobscured AGN, i.e a power-law X-ray spectrum with photon index 
T = 1.9. 



Q.M = 0.3, J1a = 0.7). Therefore, by selecting galaxies at z « 1, 
it is more likely to identify systems that are physically similar to the 
IRX population at ~ 2. In this respect it is interesting that in the 
local Universe there are no known analogues of the IRX sources, 
which can be used to study in detail the nature of this population 
(e.g. |Dey et al.|2008] >. 

Throughout this paper we adopt Ho = 70 km s _1 Mpc -1 , 
fi M = 0.3 and Q A = 0.7. 



2 DATA 

The UV to far-IR data used in this study are from the Universidad 
Complutense de Madrid (UCM) Cosmological Surveys Databas^] 
which combines all the photometric and spectroscopic observations 
available in prime cosmological survey fields. Details on the data 
reduction, band merging and photometry are presented by |Perez^| 
Gonzalez et al. ( 2005 , 2008 ). The fields of choice for this paper are 
the AEGIS and GOODS-North. In both fields the IRAC 3.6 
selected catalogues are used, which reach limiting magnitudes (AB 
system) of [3.6]=25.5 and 24.5 mag respectively. A total of 80237 
and 9061 sources are detected in the AEGIS and GOODS-North 
respectively to the limits above. 

The multiwavelength data in AEGIS include (i) photometry in 
the GALEX far- and near-ultraviolet wavebands (FUV and NUV 
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Figure 2. Distribution of X-ray-to-mid-IR luminosity ratio, 
Lx(2 — 10keV)/vL I/ (5.8/xm). The (black) unshaded histogram is 
for AGN from Lutz et al. (2004). This sample was selected to include 
only AGN that dominate in the mid-IR. The vertical dashed lines show the 
dispersion of the X-ray-to-mid-IR luminosity ratio for the same sample. 
The low-luminosity AGN of Terashima et al. (2002) are shown with the 
blue hatched histogram. The mean Lx(2 — 10 heV)/uL v (5.8/xm) of 
IRX sources is plotted with the (red) arrow. This estimate is for X-ray 
undetected IRX sources in the CDF-North (Georgantopoulos et al. 2008). 
The mean X-ray luminosity of this population is determined via stacking 
analysis and vL v (5.8/jm) is the average mid-IR luminosity of the sample. 
Both luminosities are estimated assuming 2 = 2 for the IRX sources. The 
5.8/tm luminosity of the Terashima et al. (2002) sources is determined us- 
ing the IRAS 12 fim flux density and adopting an average flux density ratio 
/12 nm//5.8 Mm = 0-9, estimated from the subsample of these sources 
with 6fim flux density measurements from the literature. The 2-10keV 
X-ray luminosities of the Terashima et al. (2002) AGN are corrected for 
obscuration, except for a total of 6 sources in that sample, which show 
a strong FeKa 6.4 keV line in their X-ray spectra with equivalent width 
EW > 900 eV and are therefore Compton Thick candidates. 



respectively; Davis et al. 2007), (ii) BRI observations from the 
Canada-France-Hawaii Telescope CFHT12K mosaic camera ( |Coil| 
|et al.|2004| >, (iii) ugriz imaging from the T0003 data release of the 
Canada-France-Hawaii Telescope Legacy Survey (CFHTLS ; |Ilbert| 
|et al.||2006|>, (iv) Spitze r IRAC observations at 3.6, 4.5, 5.8 and 
8.9 fim |Barmby et al.||200"8] > obtained as part of the Guaranteed 
Time Observing (GTO) program and (v) Spitzer MIPS photometry 
at 24 and 70 fim by coadding GTO data with the 2nd data release of 
the Far-Infrared Deep Extragalactic Legacy Survey (FIDEL). Spec- 
troscopy in the AEGIS is primarily from the DEEP2 survey ( |Davis| 
et al. 2003 , 2007 1 and a number of smaller spectroscopic programs 
that targeted the original Groth Strip ( Weiner et al. 2005 and refer- 
ences therein). 

The photometry in the GOODS-North includes (i) GALEX 
FUV and NUV data, (ii) deep UBVRIz ground based imaging 
presented by Capak et al. (2004), (iii) Spitzer IRAC (3.6, 4.5, 5.8 
and 8.9 /im) and MIPS 24 /im observations obtained as part of the 
GOODS program, MIPS 70 fim data from FIDEL. Optical spec- 
troscopy in the GOODS-North is available from either programmes 
that specifically target the X-ray population in this field (e.g. Barger 
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|et al.|2003||2005||Cowie et al.|2003| > or the Keck Treasury Redshift 
Survey (TKRS; |Wirth et al.|2004| > 

The X-ray observations of the Extended Groth Strip (AEGIS- 
X) consist of 8 ACIS-I (Advanced CCD Imaging Spectrometer) 
pointings, each with a total integration time of about 200 ks split 
in at least 3 shorter exposures obtained at different epochs. The 
GOODS-North lies in central and most sensitive region of the 
Chandra Deep Field North (CDF-N) survey, which consists of 20 
individual ACIS-I observations, which sum up to a total exposure 
time of 2Ms. The CDF-N and AEGIS-X data reduction, source 
detection and flux estimation have been carri ed out in a homoge- 
neous way using the methodology described by |Laird et al.| {2009 ^ 
Sources with Poisson probabilities < 4 x 10 ~ a are detected in 
4 energy bands, full (0.5-7 keV), soft (0.5-2 keV), hard (2-7 keV) 
and ultra-hard (5-7 keV). The count rates in the above energy in- 
tervals are converted to fluxes in the standard bands 0.5-10, 0.5- 
2, 2-10 and 5-10keV, respectively, assuming a power-law X-ray 
spectrum with T = 1.4 absorbed by the Galactic hydrogen col- 
umn density appropriate for each field. The 2-10keV fluxes where 
the sensitivity curves drop at half their maxima are 2 x 10 -14 and 
2 x 10~ 16 erg s -1 cm" 2 for the AEGIS-X and CDF-N surveys re- 
spectively. The identification of the X-ray sample with the IRAC 
3.6 fim sources has been carried out using the The Likelihood Ra- 
tio method (LR; [Sutherland & Saundersp992l |Laird et al.|2009l 
|Georgakakis et a ,|2009^ In the AEGIS, of the 918 X-ray sources 
that overlap with the Spitzer IRAC survey region, 867 (93 per cent) 
have an IRAC counterpart with likelihood LR > 0.5 and an es- 
timated contamination rate of 1.5 per cent. In the CDF-North 335 
X-ray sources overlap with the GOODS-North IRAC data, of which 
301 (90 per cent) have counterparts with LR > 0.5. The expected 
contamination rate is about 1 .4 per cent. 



3 SAMPLE SELECTION 

Lower redshift analogues of the z ~ 2 IRX population are iden- 
tified as follows. Firstly, IRAC 3.6 \rm sources with spectroscopic 
redshifts only are selected. When a spectroscopic quality flag is 
available we select sources with redshift determinations secure at 
the > 90% confidence level (i.e. quality flag Q > 3 for the DEEP2 
and the TKRS surveys). This reduces the sample to 1784 sources 
in the GOODS-North and 6488 the AEGIS. This is to avoid un- 
certainties in the determination of the rest-frame SEDs associated 
with insecure spectroscopic redshift determinations or photomet- 
ric redshift errors. The rest-frame SEDs of the spectroscopically 
identified IRAC 3.6 /im sources are then redshifted to z = 2, the 
mean redshift of the IRX population. We select only those sources 
which at z — 2 have observed SEDs with /24(im//ij > 900 and 
R— [3.6] > 3.7 mag (AB system), i.e. similar to the criteria used by 
Fioreetal. |(2 008 2009) to define their IRX sample. |G"e"organtopou^| 



los et al 



1 2008 1 have shown that the colour cut R— [3.6] > 3.7 mag 



adopted here is similar to the R — K > 4.5 mag used by |Fiore et al.| 
( |2008| |2009) . The method for selecting lower redshift analogues 
of the z w 2 IRX population is demonstrated by the dotted lines 
in Figure [5] which plots the flux density ratio /24 M m//ii against 
R - [3.6] colour. A total of 19 sources (14 in the AEGIS and 5 in 
the GOODS-North) are selected. The sample is presented in Table 
[T] The observed flux densities of the sources are shown in Table [2] 
for the AEGIS and Table[3]for the GOODS-North galaxies. 



data available at http://astro.imperial.ac.uk/research/xray/ 



For the estimation of the 2 = 2 SED we interpolate linearly 
between observed data points. It is recognised that the mid-IR is 
a complex wavelength regime with PAHs and/or the silicate ab- 
sorption/emission feature at 9.7 )im having a strong impact on the 
observed broad-band flux densities. The linear interpolation ap- 
proach is ignoring these subtleties but is preferred over more com- 
plex template fitting schemes as it is model independent. It is con- 
firmed nevertheless, that the sources presented in Table [T] satisfy 
the IRX criteria at z — 2 if the SED fits described in section [4] are 
used to interpolate between data points. The Spitzer MIPS 24 )im 
waveband and the i?-band correspond to rest-frame wavelengths of 
about 8 fim and 0.22 /im at z — 2 respectively. In order to avoid 
uncertain extrapolations of the SED we only consider sources with 
flux density detections, not upper limits, at rest-frame wavelengths 
shorter than 0.22 fim and longer than 8 /im. In practice this means 
that for the typical redshift of the sources in the sample, they should 
be detected at 24 or longer wavelengths and the U -band or 
shorter wavelengths. Three sources in the AEGIS are not detected 
at 24 /im and therefore the observed 24 (im flux density at z — 2 
is estimated by interpolating the SED between 8 and 70 /im. 

We caution that the sample of IRX sources selected in this 
paper is by no means complete. For example, the requirement for 
secure spectroscopic redshifts translates to an optical magnitude 
limit of R « 24 mag. It can be shown however, that this does 
not introduce a bias to the sample, i.e. in favour of optically bright 
sources. We calculate that the median observed i?-band flux den- 
sity of our sources if they were at z = 2 would be ss 0.04 /iJy. 
This is estimated after scaling with the appropriate k-corrections 
(i.e. 1 + 2 factor) and luminosity distances the B-band flux den- 
sities of our sources, which roughly correspond to the i?-band at 
2 = 2. The number above should be compared with the median R- 
band flux density of « 0.07 /uJy for sources in the AEGIS and 
the GOODS-North fields with observed /a4 M m//fl > 900 and 
R — [3.6] > 3.7 mag (AB system). If anything our sample is sen- 
sitive to somewhat optically fainter sources compared to the 2 « 2 
IRX population in the AEGIS and the GOODS-North fields. Ad- 
ditionally, the median 5.8/^m luminosity of our IRX sample, is 
^L„(5.8/im) w 1.6 x 10 44 ergs" 1 , about 2 times fainter than the 
median 5.8/im luminosity of the IRX sources in the CDF-South 
(S24 > 40/i Jy; Fiore et al. 2008) and 1 dex fainter than those in 
the COSMOS (S 24 > 500^Jy; Fiore et al. 2009). We therefore 
conclude that our z « 1 sample, although not complete, is similar 
to the 2 ~ 2 IRX sources found in the CDF-South, for which Fiore 
et al. (2008) estimate a Compton Thick rate of 80 ± 15 per cent and 
an intrinsic luminosity of Lx{2 — lOkeV) > 10 43 ergs _1 . 



4 SED FITTING 

The observed optical to mid-IR Spectral Energy Distribution (SED) 
of the sample sources are modeled following the methods fully 
described in Rowan-Robin son et al.| ( 2005 , 2008 ). In brief the 
FUV -band to 4.5 (im photometric data are fit using a library of 
8 templates described by Babbedge et al. (2004), 6 galaxies (E, 
Sab, Sbc, Scd, Sdm and sb) and 2 AGN. At longer wavelengths 
(5.8 — 24 /j,m) any dust may significantly contribute or even dom- 
inate the observed emission. Before fitting models to these wave- 
lengths the stellar contribution is subtracted from the photometric 
data by extrapolating the best-fit galaxy template from the previous 
step. The residuals are then fit with a mixture of four templates: cir- 
rus (Efstathiou & Rowan-Robinson 2003 1, AGN dust tori (Rowan- 
|Robinson|[T995l pfstathiou & Rowan-Robinson|[T995l >, M 82 and 
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1678 


o.u 
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1 1 7 
1.1/ 


irx-6 


14 1 f\ s 1 on 

14- lO JJ .UU 


_US? ?^ 1 7 £8 
-f-JZ Zj 1 / .Do 


1 (Y)AA 


1 1 (XX 
1 LVD 


S 8 

J.O 
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c 




1 79 
1 . /Z 




14 17 9^ 10 


-US? 77 S7 
tjz z, / ny.j / 


U.7UZM- 


1 8Q1 


D. 1 


1 1 41 








1 71 
l . / 1 


irx-8 


14 17 41.87 


+52 28 23.42 


1.1482 


1753 


5.3 


12.34 


12.26 


E 


_ 


2.64 


irx-9 


14 17 45.94 


+52 30 32.57 


0.9853 


1417 


5.2 


11.62 


11.12 


E 




1.92 


irx-10 


14 18 21.86 


+52 38 42.08 


0.7197 


1365 


5.7 


11.23 




E 




1.53 


irx-11 


14 19 38.59 


+52 55 53.69 


0.8976 


1661 


4.6 


11.31 




Sab 




1.61 


irx-12 


14 20 38.87 


+53 08 16.21 


0.6873 


1531 


5.7 


11.07 


10.77 


E 




1.37 


irx-13 


14 21 16.46 


+53 05 27.55 


0.7433 


1282 


6.0 


11.07 




E 




1.37 


irx-14 


14 22 38.98 


+53 24 14.80 


1.2873 


2984 


4.6 


12.49 




E 




2.79 


irx-15 


12 36 29.13 


+62 10 45.95 


1.0130 


1055 


5.6 


11.91 




E 




1.89 


irx-16 


12 36 56.01 


+62 08 07.86 


0.7920 


1136 


5.7 


11.78 




E 




2.08 


irx-17 


12 36 46.61 


+62 10 48.44 


0.9399 


912 


5.1 


11.41 




Sab 


0.35 


1.39 


irx-18 


12 37 16.65 


+62 17 33.56 


1.1460 


1112 


4.5 


11.48 


11.23 


Sab 


0.2 


1.78 


irx-19 


12 37 48.67 


+62 13 12.84 


0.9110 


1641 


5.5 




11.75 


QSO 


1.50 





The columns are: (1): Source identification; (2): Right Ascension (J2000) of the IRAC 3.6 mum source; (3) Declination (J2000) of the IRAC 3.6 mum 
source; (4) spectroscopic redshift; (5) observed /24pm//i? flux ratio if the source were at z = 2; (6) observed R — [3.6] colour (AB system) if the source 
were at z = 2; (7) log of the IR luminosity (3 — lOOO^tm) of the starburst component, estimated from the SED spectral fits. A dash indicates that this 
component is not required to fit the data; (8) log of the IR luminosity (3 — lOOO^tm) of the AGN torus component, estimated from the SED spectral fits. A 
dash indicates that this component is not required to fit the data; (9) best-fit optical template in the optical; (10) dust extinction; (11) log of the star-formation 
rate estimated from the SED fits. A dash indicates no dust reddening. 



Table 2. Multiwavelength SEDs of the AEGIS IRX sources 



ID 


u 


B 


3 


R 


r 






I 










i 










z 








3.6^m 




4.5/J 


m 




5.8/j 


m 


irx-1 


-0.92 ± 0.12 


-0.81 ± 0.10 


-0.33 ± 0.04 


0.45 ± 0.01 


0.31 ± 0.01 





93 


± 





01 





80 


± 





01 





97 


± 





01 


1 


63 


± 


(1 


01 


1 


43 


± 





01 


1 


35 


± 


0.(1 


irx-2 


0.10 ± 0.03 


0.04 ± 0.01 


0.23 ± 0.01 


0.43 ± 0.01 


0.48 ± 0.01 





76 


± 





01 





82 


± 





01 





95 


± 





01 


1 


88 


± 





01 


2 


06 


± 





01 


2 


22 


± 


0.0 


irx-3 


-0.90 ± 0.12 


-0.70 ± 0.10 


-0.48 ± 0.04 


0.27 ± 0.01 


0.08 ± 0.02 





74 


± 





01 





56 


± 





01 





71 


± 





01 


1 


80 


± 





01 


1 


64 


± 





01 


1 


61 


± 


0.0 


irx-4 


-0.82 ± 0.22 


-0.46 ± 0.10 


-0.21 ± 0.03 


0.41 ± 0.01 


0.32 ± 0.01 





98 


± 





01 





85 


± 





01 


1 


09 


± 





01 


2 


09 


± 





01 


1 


95 


± 





01 


1 


83 


± 


0.0 


irx-5 


-0.86 ± 0.14 


-0.36 ± 0.04 


0.03 ± 0.02 


0.75 ± 0.01 


0.64 ± 0.01 


1 


16 


± 





01 


1 


09 


± 





01 


1 


23 


± 





01 


1 


90 


± 





01 


1 


85 


± 





01 


1 


89 


± 


0.0 


irx-6 


-0.83 ± 0.10 


-0.83 ± 0.16 


-0.41 ± 0.04 


-0.04 ± 0.02 


-0.01 ± 0.02 





54 


± 





01 





40 


± 





01 





65 


± 





01 


1 


72 


± 





01 


1 


64 


± 





01 


1 


44 


± 


0.0 


irx-7 


-0.84 ± 0.10 


-0.98 ± 0.20 


-0.38 ± 0.03 


0.15 ± 0.03 


0.00 ± 0.02 





65 


± 





01 





53 


± 





01 





74 


± 





01 


1 


75 


± 





04 


1 


61 


± 





04 


1 


48 


± 


0.0 


irx-8 


0.26 ± 0.02 


0.11 ± 0.01 


0.41 ± 0.01 


0.93 ± 0.01 


0.84 ± 0.01 


1 


39 


± 





01 


1 


27 


± 





01 


1 


49 


± 





01 


2 


69 


± 





01 


2 


89 


± 





01 


3 


05 


± 


0.0 


irx-9 


-0.71 ± 0.14 


-0.35 ± 0.04 


-0.25 ± 0.03 


0.38 ± 0.02 


0.32 ± 0.01 





98 


± 





01 





85 


± 





01 


1 


13 


± 





01 


1 


92 


± 





01 


1 


92 


± 





01 


2 


01 


± 


0.0 


irx-10 


-1.05 ± 0.16 


-0.60 ± 0.08 


-0.44 ± 0.04 


0.23 ± 0.02 


0.02 ± 0.02 





59 


± 





01 





45 


± 





01 





63 


± 





01 


1 


68 


± 





01 


1 


51 


± 





01 


1 


5 5 


± 


0.0 


irx-11 


-0.86 ± 0.13 


-0.70 ± 0.11 


-0.33 ± 0.03 


0.12 ± 0.01 


0.07 ± 0.02 





56 


± 





01 





45 


± 





01 





63 


± 





01 


1 


31 


± 





04 


1 


10 


± 





04 





94 


± 


0.0 


irx-12 


-0.58 ± 0.06 


-0.46 ± 0.08 


-0.06 ± 0.02 


0.62 ± 0.01 


0.49 ± 0.01 


1 


07 


± 





01 





92 


± 





01 


1 


09 


± 





01 


2 


03 


± 





01 


2 


03 


± 





01 


2 


09 


± 


0.0 


irx-13 


-1.39 ± 0.45 


-0.58 ± 0.16 


-0.55 ± 0.04 


0.38 ± 0.01 


0.06 ± 0.02 





86 


± 





01 





57 


± 





01 





76 


± 





01 


1 


71 


± 





01 


1 


54 


± 





01 


1 


51 


± 


0.0 


irx-14 




-0.43 ± 0.06 




0.08 ± 0.02 







51 


± 





04 






















1 


75 


± 





04 


1 


80 


± 





04 


1 


65 


± 


0.0 



None of the sources are detected in the GALEX FUV band. Source irx-8 has a GALEX NUV log flux density of 0.10 ±0.03 (/xjy units). The columns are: (1): Source identification; (2): log of CFHTLS 
u-band flux density in units of /i Jy; (3): log of DEEP2 S-band flux density in units of //Jy; (4): log of CFHTLS p-band flux density in units of /j Jy; (5): log of DLEP2 _R-band flux density in units of //Jy; 
(6): log of CFHTLS r-band flux density in units of [i Jy ; (7): log of DEEP2 /-band flux density in units of fi Jy ; (8): log of CFHTLS i-band flux density in units of fj, Jy; (9): log of CFHTLS z-band flux density 
in units of fj, Jy; (10): log of IRAC 3.6//m-band flux density in units of (j, Jy; (11): log of IRAC 4.5/J.m-band flux density in units of /iJy; (12): log of IRAC 5.8/j.m-band flux density in units of fi Jy; (13): 
log of IRAC 8.9/im-band flux density in units of fJ,Jy; (14): log of MIPS 24//m-band flux density in units of fiJy ; (15): log of MIPS 70/j.m-band flux density in units of fi Jy. 



Table 3. Multiwavelength SEDs of the CDF-North IRX sources 



ID 


U 


B 


R 


I 


z 


3.6^tm 


4.5^m 


5.8//m 


8.0Mm 


24nm 


70jum 


irx-15 
irx-16 
irx-17 
irx-18 
irx-19 


-0.69 ± 0.05 
-0.70 ± 0.04 
-0.83 ± 0.05 
-0.44 ± 0.02 
-0.94 ± 0.08 


-0.52 ± 0.04 
-0.30 ± 0.03 
-0.75 ± 0.05 
-0.33 ± 0.04 
-0.80 ± 0.05 


0.11 ± 0.01 
0.57 ± 0.01 
-0.04 ± 0.01 
0.22 ± 0.01 
0.13 ± 0.01 


0.59 ± 0.01 
1.00 ± 0.01 
0.39 ± 0.01 
0.57 ± 0.01 
0.64 ± 0.01 


0.78 ± 0.01 
1.17 ± 0.01 
0.54 ± 0.01 
0.71 ± 0.01 
0.83 ± 0.01 


1.91 ± 0.01 
2.00 ± 0.01 
1.49 ± 0.01 
1.73 ± 0.01 
1.56 ± 0.01 


1.87 ± 0.01 

1.85 ± 0.01 
1.37 ± 0.01 

1.86 ± 0.01 
1.65 ± 0.01 


1.84 ± 0.01 
1.88 ± 0.01 
1.39 ± 0.04 
2.07 ± 0.01 
1.91 ± 0.01 


1.80 ± 0.02 
1.89 ± 0.02 
1.34 ± 0.04 
2.31 ± 0.01 
2.12 ± 0.01 


2.73 ± 0.01 
2.85 ± 0.01 
2.47 ± 0.01 
2.99 ± 0.01 
2.53 ± 0.01 


3.24 ± 0.07 
4.09 ± 0.02 
3.39 ± 0.05 
3.41 ± 0.06 



None of the sources are detected in the GALEX NUV and FUV bands. The columns are: (1): Source identification; (2): log of £/-band flux density in units of //Jy; (3): log of S-band flux density in units of 
fi Jy; (4): log of DEHP2 H-band flux density in units of /z Jy; (5): log of DELP2 /-band flux density in units of /i Jy; (6): log of CFHTLS z-band flux density in units of //Jy; (7): log of IRAC 3.6//m-band 
flux density in units of //Jy; (8): log of IRAC 4.5//m-band flux density in units of //Jy; (9): log of IRAC 5.8//m-band flux density in units of //Jy; (10): log of IRAC 8.9//m-band flux density in units of 
//Jy; (11): log of MIPS 24//nvband flux density in units of //Jy; (12): log of MIPS 70//m-band flux density in units of //Jy; 
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Figure 4. Examples of template fits to the SEDs of the IRX selected galaxies in the AEGIS. From bottom to top are: irx-8 (X-ray source), irx-9 (X-ray source), 
irx-10, irx-11, irx-12 (X-ray source), irx-13, irx-14 and irx-7. The dots are the observed UV-to-far-IR photometry and the continuous lines are the best-fit 
models. Different sources are offset by an arbitrary constant for clarity. The dashed curves correspond to the hot dust torus component, which is needed for all 
X-ray detected sources with Lx (2 — 10 keV) > 10 43 erg s _1 . The dotted lines at infrared wavelengths are the starburst components. 



Arp220 starbursts (Efstat hlou et al.| |2000). The modeling above 
provides both information on the dominant emission mechanism 
in the optical and the infrared (AGN vs star-formation) and an es- 
timate of the total infrared luminosity, Ltot, in the wavelength 
range 3 — 1000/im. As discussed by Rowan-Robi nson et al.H2005) 
Ltot is expected to be accurate within a factor of two. In this exer- 
cise the redshift is fixed to the spectroscopically determined value. 

The results are presented in Table [T] Examples of template 
fits to SED of AEGIS and CDF-North sources are shown in Fig- 
ures [4] and [5] Our analysis shows that dust associated with star- 
formation (M82 starburst template) either dominates or substan- 
tially contributes to the flux at the mid/far-infrared wavelengths of 
most sources. In this respect, it is not surprising that 12/14 galax- 
ies in AEGIS and 4/5 sources in the CDF-N are detected at this 
longer wavelength. The estimated star-formation rates are in the 
range 10 — 600 Mq /yr. The best-fit starburst template is that of 
M 82, in agreement with results from Buss mann et al.|j2009a[ > who 
showed that the far-IR SEDs of IRX sources are better described 



by warm-dust templates (e.g. similar to M82 or Mrk231) and are 
inconsistent with cold dust templates (e.g. Arp 220). 

An additional QSO torus component is required to fit the mid- 
IR part of the SED of some sources in Table [T] As will be shown 
in the next section this component is needed for only those sources 
which are detected at X-ray wavelengths and which have intrinsic 
X-ray luminosities, L x (2 - lOkeV) > 10 43 ergs" 1 . The SED fit- 
ting approach can therefore be used to identify AGN with intrinsic 
luminosity > 10 43 ergs -1 . Below this limit, star-formation in the 
host galaxy may swamp any AGN signatures in the mid-IR. 



5 X-RAY PROPERTIES 

5.1 X-ray source counterparts 

The cross-correlation between the X-ray and IRAC 3.6/im cata- 
logues using the Likelihood Ratio method shows that 5/16 and 4/5 
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Figure 5. Same as in figure[4]for IRX selected galaxies in the CDF-North. From bottom to top are: irx-15, irx-16, irx-17, irx-18, irx-19. The curves and dots 
are the same as in Figure[4] 



of the sample sources in AEGIS and the CDF-North respectively, 
have X-ray counterparts with Poisson significance < 4 x 10 in 
at least one of the 4 energy bands used for detection. Lower sig- 
nificance X-ray counterparts are also searched for in the hard band 
(2-7 keV), which provides a reasonable compromise between sensi- 
tivity and obscuration effects, thereby minimising biases when ex- 
trapolating fluxes to the rest-frame 2-10keV energy range. The X- 
ray counts within the 70 per cent Encircled Energy Fraction (EEF) 
radius at the position of each source in Table 1 are summed up 
and a local value for the background is estimated using an annulus 
centered on the source with inner radius 1.5 times the 90 per cent 
EEF and width of 50arcsec (100 pixels). For the estimation of the 
background X-ray sources with Poisson significance < 4 x 10~ 6 
are removed by excluding pixels within the 95 per cent EEF ra- 
dius of each source. The probability that the observed counts are 
a random fluctuation of the background is estimated using Poisson 
statistics. We consider as detections sources with Poisson signifi- 
cance < 3 x 10~ 3 . This exercise yields 2 lower significance X-ray 
detections, 1 in the AEGIS and 1 in the CDF-N. The fluxes of these 
sources are estimated using a Bayesian methodology similar to that 



described in Laird et al. ( 2009 ). For a source with observed number 
of total counts N (source and background) within the 70 per cent 
EEF radius and a local background value B the probability of flux 

fx is 



P(fx,N) = 



Nl 



Ah 



(l) 



where N = S + B and S = fx x t exp x C x r\. In the last equation 
t exp is the exposure time at a particular position after accounting 
for instrumental effects, C is the conversion factor from flux to 
count rate and rj is the encircled energy fraction, i.e. 0.7 in our case. 
The term ir(fx) in equation[T|accounts for the Eddington bias and 
assumes that the differential X-ray source counts follow a broken 
power-law for the log iV — log S with faint and bright-end slopes of 
-1.5 and -2.5 respectively and a break flux of 10 _14 ergs _1 cm -2 
(Georgakakis et al. 2008). For the source flux we adopt the mode 
of the distribution P(fx, N). 

The same methodology is also adopted to estimate upper lim- 
its for the fluxes of the sources that are not detected at X-ray wave- 
lengths. The confidence interval CL is given by the integral 
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Table 4. X-ray properties of the sample 



ID XID counts bkg N H R V c - stat dof L x (2 - 10 keV) 

! cm- 2 ) (xlO 43 ergs" 1 ) 



irx-1 

irx-2 

irx-3 

irx-4 i 

irx-5 

irx-6 

irx-7 

irx-8 

irx-9 

frx-10 

irx-1 1 

irx-1 2 

irx-1 3 

irx-1 4 

irx-1 5 
irx-1 6 
irx-17 
irx-1 8 
irx-19 



egs_0263 



egs_0311 



egs_0416 
egs_0424 



egs_0917 



hdfh_098 
hdfn_328 
hdfn_466 
hdfn_023 



16 

127 

4 

112 
14 

2 
820 
254 



6 

163 
19 

2 

174 
129 

27 
2240 
260 



7.3 
0.6 
0.2 
0.5 
5.9 
10.4 
0.2 
1.0 

I. 2 
0.1 
4.1 
47.0 

II. 8 
0.2 

28.0 
96.7 
20.5 
74.2 
210.1 



44.C 



7.9 



h26.0 
-12.1 



-2.!) 
-2.5 



0.6 



+1.2 

0.6 



2.1 

4.6 



+0.6 
0.8 
1.1 
1.1 



i0.3+J; 5 5 



9 5+ 57 

1 3+ 4 1 

< 2.8 

o n+0.8 
z -°-0.6 



2.3 



+8.0 

2.9 



03 +7 h 2 I 1 
u - uo -0.03 z -0.2 

0.5^-7 1.9(f) 



55.4 



-117.5 
-38.2 



6.4 



5.6 
0.2 



1.9(f) 



1.9(f) 



,+0.7 



2.6 



42.1 



40.5 



48.2 
44.9 



46 



46 



45 
46 



< 0.7 
36.4 

< 0.07 
0.1 
3.2 

< 0.6 

< 0.9 
27.9 











< 0.1 










< 0.4 


+ 10.7 
'-2.6 


1.9(f) 


43.4 


46 


2.4 










< 0.4 










< 1.2 


+22.8 
-6.4 


1.9(f) 


41.4 


39 


0.6 


1(f) 


1.9(f) 


50.3 


47 


0.05 


1(f) 
+22.5 
'-2.5 
,+18.4 
'-0.2 


1.9(f) 
2 0+ ' 7 
1.9(f) 


43.1 
51.4 
103.5 


49 
45 
117 


0.046 
9.3 
0.7 



The columns are: (1): Source identification; (2): X-ray identification number as in Laird et al (2009) for AEGIS and Georgakakis et al. (2008) for the CDF- 
North. Sources without an X-ray counterpart or low significance (Poisson probability > 4 X 10~ 6 ) X-ray sources do not have an X-ray identification number 
(3): total counts in the energy range 0.5-10keV. For sources with insufficient counts to perform X-ray spectral analysis the total extracted counts in the 2-8 keV 
band are listed. (4): background counts in the range 0.5-10keV. For sources with insufficient counts to perform X-ray spectral analysis the background counts 
in the 2-8 keV band are listed. (5): best-fit column density of the cold absorber (WABS component of XSPEC model). The uncertainties correspond to the 90 
per cent confidence level; (6): reflection fraction of the PEXMON model of XSPEC. The errors correspond to the 90 per cent confidence level; (7): power-law 
photon index, T, of the intrinsic AGN power-law spectrum. The "(f)" indicates that this parameter was kept fixed in the fit. (10): C-stat of the best-fit solution; 
(9): degrees of freedom; (11) obscuration corrected luminosity in the 2-10 keV range and in units of 10 43 erg s — . In the case of upper limits or sources with 
insufficient counts to perform X-ray spectral analysis, no correction has been applied for intrinsic obscuration. In those sources fluxes and luminosities in the 
2-10keV band are estimated by adopting a power-law spectral index with T = 1.4, i.e. assume that the sources have X-ray spectra similar to that of the X-ray 
background. The statistical uncertainty in the luminosity is less than 10% for X-ray detected sources (Poisson probability < 4 X 10~ 6 ) and about a factor of 
2 for low significance ones. 



[ fX U P(fx,N)df x =CL. (2) 

Jfx.L 

For a faint-end slope of -1.5 the integral diverges at the 
faintest fluxes. To avoid this problem we adopt as the lower limit 
of the integration fx,L = 1CP 18 ergs -1 cm~ 2 . Equationj^jis then 
solved numerically to estimate the flux upper limit fx.u for the 
confidence interval of 99.97 per cent. The results are presented in 
Table 1. 



5.2 X-ray spectral analysis 

The X-ray spectra of the sample sources were extracted using the 
ACIS extract (AE; v. 2008-03-04; |Broos et al.|2002) > IDL package, 
which is designed to deal with large number of sources on multi- 
ple observations. The AE package User's Guide is available online 
at http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide.html. 
The procedures used in AE are described in Townsley et al.|p003) 
and |Getman et al.] < [2005| >. The X-ray counts of each source were 
extracted from individual pointings using the 95 per cent encir- 
cled energy radius (1.5 keV) at the position of the source. The 
composite source spectrum was then constructed by summing the 
counts from each observation. For the background estimation, point 
sources were first masked out using a circular region with size 1.1 
larger than the 99 per cent EEF at 1.5 keV. The background spec- 



tra were then extracted from the masked event files of individual 
observations using an annulus centered on the source with vari- 
able size that enclosed at least 150 background counts. The spectra 
from separate pointings were scaled by the ratio of the exposure 
time in the source and background regions and then summed up 
into the composite background spectrum for a particular source. 
The Response Matrix Files (RMF) and Area Response Files (ARF) 
for individual observations were constructed using the CIAO tools 
MKACISRMF and MKARF respectively. These were then combined 
into composite RMFs and ARFs by the FTOOL tasks ADDRMF and 
ADDARF respectively, using the exposure time of each observation 
as weight. Source spectra were grouped to have at least one count 
per bin. 

For the spectral fits we use the XSPEC vl2. We adopt a realis- 
tic AGN spectral model to fit the data, motivated by current obser- 
vations and ideas about the structure of the regions close to SBHs. 
The model consists of power-law absorbed by cold material and a 
reflection component as described below. 

The main spectral component of AGN is the power law with 
an average photon index F — 1.9 (Nandra & Pounds 19941 and a 
cutoff at high energies (~ 100 keV; Gondek et al. 19961. The dis- 
persion around the mean F is estimated to be 0.2-0.3 ( jGondek et af] 
1996). The reflection component represents direct AGN emission 
that intersects optically thick material, such as the accretion disk or 
the torus, and is Compton backscattered into the line of sight. We 
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Figure 3. /24/jm//jj against R — [3.6] estimated by redshifting to z = 2 
the rest-frame SEDs of 3.6 fim sources in the AEGIS and the CDF-North 
with spectroscopic redshifts (blue dots). Large red circles are the subset of 
those sources which at z = 2 have /m^ih/ }r > 900 and R — [3.6] > 
3.7 mag (dotted horizontal and vertical lines respectively), i.e. similar to 
the IRX source selection criteria (Fiore et al. 2008, 2009; Georgantopoulos 
2009). 



parametrise this component using the PEXMON model of XSPEC 
(Nandra et al. 2007b I, which includes an analytic formulation for 
the strength of the narrow iron line at 6.4 keV, which is expected 
to be prominent in the reflection spectrum. The direct emission of 
AGN is also modified by cold material. We adopt the photoelectric 
absorption cross sections of Morrison & McCammon (1983) for 
solar metallicity to represent any cold absorbing medium, such as 
a torus. Ideally, our modeling should also include absorption by 
ionised material. However, given the relatively small number of 
counts in the typical X-ray spectra of deep field sources, we choose 
not to include this component in the analysis. In any case, warm 
absorption is predominantly affecting the AGN spectra at soft en- 



ergies < 1 keV (e.g. George et al. 1998 1. At least part of the warm 
absorber spectral features are therefore, expected to be redshifted 
out of the Chandra observable energy range for sources at mod- 
erate and high redshifts. We also perform the spectral modelling 
for source photons > 2 keV (observer's frame) to ensure that soft 
emission associated with a warm absorber does not modify the re- 
sults and conclusions. In summary the XSPEC model we are using 
consists of a power-law absorbed by cold gas at the rest-frame of 
the source and Compton reflection. These components are also ab- 
sorbed by cold gas in our Galaxy using the appropriate Galactic HI 
column density for the AEGIS (Nh = 1-4 x 10 20 cm~ 2 ) and the 
CDF-North (N H = 1.3 x 10 20 cm' 2 ). In XSPEC terminology, our 
model is wabs X (zwabs X po + pexmon), where wabs and zwabs 
is the cold absorption from our Galaxy and at the rest frame of the 
sources respectively, and po is the power-law. 

The parameters of the model are the photon index and the high 
energy cutoff (pexmon only) of the direct AGN emission, the col- 
umn density of the cold absorber, Nh, the element abundances, the 
inclination i of the slab that gives the reflection and the reflection 
fraction R, defined as the strength of the reflection relative to that 
expected from a slab subtending a solid angle of 27r. Clearly, these 



Figure 6. X-ray 2 — 10 keV luminosity plotted against 5.8/im luminosity. 
Filled circles (black) are local AGN from the sample of Lutz et al. (2004). 
The dashed lines show the dispersion around the mean Lx (2 — 10 keV) — 
uL v {bjj,m) relation for the Lutz et al. (2004) AGN. The (blue) long- 
dashed-dotted line is the Lx — ^Is^m relation for star-forming galaxies 
adapted from (Ranalli et al. 2003; see Fiore et al. 2009). The dotted lines 
correspond to the 1 sigma rms envelope around this relation. Our sample of 
IRX sources is shown with the (red) crosses. For IRX sources that are not 
detected at X-ray wavelengths the upper limit in Lx (2 — 10 keV) is plot- 
ted. A square on top of a cross signifies IRX sources which require a torus 
component to fit their infrared SEDs. Above an intrinsic X-ray luminos- 
ity of 10 43 erg s — 1 (shown with the horizontal green line) all IRX sources 
with an X-ray counterpart in the 2-10 keV band require a torus component. 



are too many free parameters given the number of counts in our 
spectra (see Table |4). We therefore fix abundances to solar values 
and the inclination angle to 60 deg. Given that the spectra cover the 
energy range 2-10keV (rest frame 4-20 keV at z ~ 1), the energy 
cutoff of the direct AGN emission will have negligible impact on 
the results. We therefore choose to fix this parameter to lOOOkeV 
(pexmon only). These assumptions leave 3 free parameters, T, Nh 
and R. For X-ray sources in the sample with < 200 net counts in 
the 0.5-10 keV range we also fix the photon index of the power-law 
to r = 1.9, while for sources with < 50 net counts we further fix 
the reflection fraction to R — 1. 

In the case of Compton Thick AGN we expect either to see di- 
rectly the turnover in the spectrum at rest-frame energies ^15 keV 
because of the cold absorption (i.e. estimate a column density 
Nh > 2 x 10 24 cm -2 from the spectral fit) or get a reflection 
dominated spectrum with R >> 10 (Nandra et al. 2 007a} . 

Spectral analysis was performed for sources with sufficient 
number of net counts to place meaningful constraints. This require- 
ment excluded from the spectral fits sources with low significance 
X-ray counterparts, except one, irx-19. This source lies just below 
the formal detection threshold of the CDF-North source catalogue. 
The results of the spectral analysis are summarised in Table [4] 
Spectral fits for individual sources are shown in Figures [7] [8] All 
sources in the sample are consistent with moderate column densi- 
ties in the range N H w 10 22 - 5 x 10 23 cm -2 , while the best-fit 
relative reflection parameter is < 10 for all of them. For 4 sources 
however, the 90 per cent upper limit in 7? is greater than 10. For 
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one of them, irx-19, the 90 per cent upper limit for the column den- 
sity is Nh = 1.2 x 10 24 cm -2 , close to the Compton Thick limit 
(N H = 2 x 10 24 cm" 2 ). Another of the four sources, irx-18, also 
shows the FeKa 6.4 keV in the X-ray spectrum with a rest-frame 
equivalent width of EW = 0.5 ±0.2 keV. Although a strong Fe Ka 



line (EW > 1 keV; Ghisellini et al.|l994 1 is evidence for Comp- 
ton Thick AGN, the estimated tiW is also consistent with Compton 
thin obscuration. We conclude that the X-ray spectral analysis sug- 
gests that two sources in the sample show tentative evidence for 
Compton Thick obscuration. 

Also, the X-ray emission of two X-ray detected sources (irx- 
16 and irx-17) is likely associated with starburst activity, not ac- 
cretion on a central supermassive black hole. These two sources 
have L x (2 - lOkeV) « 5 x 10 41 ergs"\ relatively soft, al- 
beit noisy, X-ray spectra and follow the correlation between X-ray 
and infrared luminosities for star-forming galaxies (Ranal lTet al.| 
|2003[ > as shown in Figure [6] Their infrared SEDs are also con- 
sistent with starburst templates and do not require an additional 
hot dust component associated with the AGN torus. Figure [6] also 
demonstrates that a QSO torus component is required to fit the mid- 
IR part of the SED of only those sources in Table [T] which are 
detected at X-ray wavelengths and which have X-ray luminosity, 
L x (2 - lOkeV) > 10 43 ergs -1 . The SED fitting method of sec- 
tionals therefore sensitive to AGN with intrinsic X-ray luminosity 
above this limit. 



6 DISCUSSION & CONCLUSIONS 

The population of IRX galaxies has attracted much attention re- 
cently as they are proposed as good candidates for Compton thick 
QSOsatz « 2 with luminosities L x (2-10 keV) > 10 43 ergs -1 . 
The low X-ray-to-mid-IR luminosity ratio and the mean hardness 
ratios of these sources are consistent with this interpretation as 
long as the bulk of the mid-IR luminosity is associated with the 
dust heated by the central engine. There are indications however, 
that this is not the case, at least not for all IRX sources. Murphy 
|et al.| ( (2009[ l combined mid-IR spectroscopy with far-IR sub-mm 
data to show that only 50 per cent of the infrared excess sources 
selected in a way similar to that described by Daddi et al. (2007), 
show evidence for obscured AGN activity. The other half of this 
population have SEDs consistent with star-formation. |Yan et aE] 
( |2007) > presented Spitzer-IRS spectra of 52 galaxies brighter than 
1 mjy at 24/xm selected in the Spitzer extragalactic First Look 
Survey (xFLS). The relative strengths of the AGN and starburst 
components in those sources have been estimated by Sajina et al. 
(2008) using in addition to the mid-IR spectroscopy, multiwave- 
length photometric observations (UV to far-IR and radio), optical 
and NIR spectroscopy. There are 21 sources in the |Yan et al.H2007) 
sample with properties similar to the IRX population at z ~ 2, 
i.e. /24nm//fl > 900 and spectroscopic redshift determination 
z > 1.5. |Sajina et al.| < |2008") found that about 76 per cent (16/25) 
of those sources are AGN dominated in the IR, while the remain- 
ing 24 per cent have starburst components that either dominate or 
contribute nearly equally to the AGN at IR wavelengths. If star- 
formation contributes to the mid-IR luminosity of IRX galaxies 
then the two key properties of this population, X-ray-to-mid-IR 
luminosity ratio and mean hardness ratio, are also consistent with 
either lower-luminosity Compton thin AGN (e.g. Figures 1, 2) or 
pure starbursts ( Donley et al. 2008}. 

In order to shed more light in the nature of IRX sources we 
select galaxies at z ~ 1 in the AEGIS and CDF-North fields with 



SEDs similar to the z ~ 2 IRX population. The advantage of this 
approach is that the selected sources have fluxes that are brighter at 
almost any wavelength compared to IRX galaxies at z ~ 2, thereby 
greatly facilitating their study. 

It is interesting that about 35 per cent of the sources in the 
AEGIS (5/14) and 80 per cent in the CDF-North (4/5) are as- 
sociated with hard (2-7 keV) X-ray detections with significance 
< 4 x 10" 6 . These fractions increase to 43 (6/14) and 100 (5/5) 
per cent in the AEGIS and CDF-N respectively, if lower signifi- 
cance sources are included. For comparison only about 4 per cent 
of red cloud galaxies in the AEGIS have X-ray counterparts (Nan- 
dra et al. 2007). This suggests a high AGN identification rate in 
samples selected using the infrared excess criteria. 

Excluding X-ray detected starburst candidates (irx-16 and irx- 
17), the IRX sources with X-ray counterparts are indeed AGN, 
as indicated by their X-ray luminosities, Lx(2 — lOkeV) > 
10 42 erg s -1 cm -2 , which are higher than what can be attributed to 
star-formation (Georgakakis et al. 2007). These sources also have 
hard X-ray spectral properties, which are however consistent with 
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There is only tentative evidence for Compton Thick obscuration 
among the X-ray detected IRX sources. For four of them the 90 per 
cent upper limit of the reflection fraction is higher than R = 10. 
Also, the X-ray spectra of two of those four sources either have 90 
per cent upper limit in the column density w 1.2 X 10 24 cm -2 or 
show the FeKa6.4keV line with an equivalent width of 0.5keV. 
These properties can be interpreted as evidence for Compton Thick 
AGN but are also consistent with Compton thin obscuration. 

More relevant to the deeply buried AGN picture are the in- 
frared excess sources in the sample that are not detected at X-ray 
wavelengths. The SED modelling shows that their mid- and far- 
IR is dominated by star-formation and that there is no indication 
for a hot dust component associated with a powerful AGN that re- 
mains undetected at X-ray wavelengths. It is noted the QSO torus 
template is required to fit the IR SEDs of all X-ray detected AGN 
with intrinsic X-ray luminosities brighter than Lx (2—10 keV) ~ 
10 43 ergs -1 . Therefore, if some of the X-ray undetected sources 
were associated with heavily obscured and powerful AGN with X- 
ray emission suppressed by the intervening dust and gas clouds, 
we would have identified them in the mid-IR as sources with a 
QSO torus component. We can therefore, place an upper limit of 
L x (2 - lOkeV) = 10 43 ergs -1 to the intrinsic AGN luminosity 
of X-ray undetected sources, if they host an active SBH. We con- 
clude that the IRX sources without X-ray counterparts in our sam- 
ple are either starbursts, lower luminosity AGN, or a combination 
of the two. 

This result has implications on the nature of the z ~ 2 IRX 
population detected in deep surveys. There is no doubt that some 
of these sources are Compton Thick QSOs. Georgantopoulo s et al.| 
(2009) for example, identified Compton Thick QSOs among X- 
ray sources in the CDF-North through X-ray spectroscopy and 
showed that some of them satisfy the infrared excess selection cri- 
teria. Our analysis however, shows that a potentially large frac- 
tion of the IRX population at z « 2 are not luminous [Lx(2 — 
lOkeV) > 10 43 ergs -1 ] Compton Thick QSOs but lower lumi- 
nosity [Lx(2 - lOkeV) < 10 43 ergs -1 ] possibly Compton thin 
AGN and/or starbursts. 

Dust enshrouded star-formation has already been proposed as 
an alternative to Compton Thick QSOs to explain the properties of 
the z w 2 IRX sources. Donley et al. (2008) critically reviewed 
different methods proposed in the literature for finding AGN at in- 
frared wavelengths, including the infrared excess selection. They 
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Figure 7. X-ray spectra and spectral fits for the IRX sources in the AEGIS. 



showed that the red SED of these sources are consistent with those 
of local pure starbursts with some moderate amounts of additional 
reddening at optical wavelengths (Ay ~ 1 mag). They also cau- 
tioned that the stacked X-ray signal of this population may be dom- 
inated by few sources. Moreover, they argued that flat mean X-ray 
spectral properties of the IRX population could be the result of high 
mass X-ray binaries, which are known to have hard X-ray spectral 
properties (power-laws with V w 1.2). Whether binary stars can 
dominate the integrated X-ray emission of local starbursts however, 
is still under debate. In any case, based on the arguments above 
|Donley et al.|(2008) suggested that as much as half of the IRX pop- 
ulation at z « 2 are starbursts with only about 20 per cent showing 
evidence for heavily obscured, possibly Compton Thick, AGN ac- 
tivity. [Pope^Ta!] {2008} analysed the mid-IR spectra of 12 IRX 
sources with f2iiim > 300 fiJy in the GOODS-North and found 



that 6 of them are dominated by star-formation at the mid-IR. They 
also showed that AGN and starburst dominated IRX sources have 
distinct 8.0/im over 4.5^im flux ratios, /s.o^m/A.s/im- Extrapolat- 
ing their results to fainter IRX sources /24 M m = 100 — 300 /iJy, 
for which mid-IR spectroscopy is not available, they found that 
about 80 per cent of them have fs.o^m/ f4.5fim consistent with 
star-formation. They also estimate the mean SED of IRX sources 
in their sample and conclude that less than about 10 per cent of 
the total infrared luminosity (8 — 1000/im) is associated with hot 
dust, possibly heated by an AGN . For the average infrared lumi- 



nosity of the Pope et al. 2008} sample, Lj R = 10 12 L , the 



fraction above translates to an upper limit in the AGN luminos- 
ity of Lir k 4 x 10 44 ergs _1 . Adopting the AGN bolomet- 
ric correction factors Lboi/Lm = 3 (Risaliti & Elvis 2004} and 
L bo i/L x (2 - lOkeV) = 35 ( jElvis et al.||1994) , we esljmate a 
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Figure 8. X-ray spectra and spectral fits for the IRX sources in the CDF-North. 



mean hard X-ray luminosity Lx (2 — 10 keV) < 3 x 10 43 erg s _1 . 
Although there are uncertainties in this calculation, there is broad 
agreement in the upper limits in Lx estimated for the Po pe et aT] 
( 2008 1 and our sample of IRX sources. 

The sample presented here has properties similar to moder- 
ately bright IRX sources at z ~ 2, i.e. 5.8 /im luminosity few 
times 10 44 ergs _1 , like those found in the CDF-South by Fiore 
et al. (2008). The simulations of Narayanan et al. (20091 suggest 
that IRX sources above 524 = 300 /xJy are dominated by gaseous 
mergers and include a large fraction of powerful AGN, whereas less 
luminous systems, like those studied here, are typically secularly 
evolving galaxies dominated by star-formation and with only weak 
AGN activity. Our results on the nature of moderately luminous 
IRX sources are therefore consistent with those simulation. It is 
likely that IRX sources more luminous than those studied here (i.e. 



COSMOS field, S 2i > 500 ^Jy, uL u (5.8fim) > lO^ergs^ 1 ; 
Fiore et al. 2009) include a higher fraction of Compton Thick 
QSOs. The results of Sajina et al. (2008) on the nature of bright 
(S24 > 900 /LtJy) IRX sources selected in the xFLS (Yan et al. 
2007) suggests that this may be the case. Bauer et al. (2009) ex- 
plored the X-ray properties of a subset of the Yan et al. (2007) 
sample and found evidence for at least mildly Compton Thick 
(Nh ss 10 24 cm -2 ) obscuration in a large fraction of the sources 
that are AGN dominated in the mid-IR. Even at the extreme lu- 
minosities of the xFLS sample however, a non-negligible fraction 
of the IRX sources (about 24 per cent) have a substantial or even 
dominant starburst component in the mid-IR. This underlines the 
importance for subtracting the contribution of star-formation to the 
mid-IR to assess the intrinsic AGN luminosity before concluding 
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whether the X-ray/mid-IR properties of these systems are consis- 
tent with Compton Thick obscuration. 

Mulla ney et al.| j2009) have recently found that the X-ray 
to total infrared luminosity ratios of X-ray AGN with Lx(2 — 
lOkeV) = 10 42 — 10 43 ergs _1 increases by about 1 dex from 
z = to z ~ 1 — 2. Although the origin of this trend is unclear, it 
can be interpreted as an enhancement of the average star-formation 
rate in lower-luminosity AGN at z ~ 2 compared to the local Uni- 
verse. In this picture, the enhanced star-formation is likely to have 
an impact on the mid-IR part of the SED and will make lower lumi- 
nosity AGN at z w 1 — 2 appear underluminous at X-rays for their 
mid-IR luminosity compared to local AGN samples. This would be 
in agreement with our interpretation of the IRX population and em- 
phasises the need to properly model the mid-IR part of the SED to 
assess the relative contribution of AGN and star-formation at this 
wavelength regime. 

Further progress in the study of IRX galaxies is expected from 
upcoming Herschel observations, which when combined with the 
existing Spitzer data will constrain the SEDs of individual sources 
at z m 2 over a sufficiently large wavelength baseline (mid to far- 
IR) to allow decomposition of the AGN and starbursts components 
through template fits, as done in this paper. This exercise will pro- 
vide estimates of the level of dusty star-formation activity and the 
intrinsic AGN luminosity of individual sources to confirm or re- 
fute claims that they are Compton Thick QSOs. The determina- 
tion of the column density of these sources through X-ray spec- 
troscopy has to wait future X-ray missions with large collecting 
areas, such IXO, which is expected to provide spectra for z w 2 
IRX sources which h ave typical X-ray fluxes fx (2 — 10 keV) 
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2009) and are cur- 



rently accessible only through stacking analysis. 
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